The strong-coupling limit of a Kondo spin coupled to a mesoscopic quantum dot: 
effective Hamiltonian in the presence of exchange correlations 
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We consider a Kondo spin that is coupled antiferromagnetically to a large chaotic quantum dot. 
Such a dot is described by the so-called universal Hamiltonian and its electrons are interacting via 
a ferromagnetic exchange interaction. We derive an effective Hamiltonian in the limit of strong 
Kondo coupling, where the screened Kondo spin effectively removes one electron from the dot. We 
find that the exchange coupling constant in this reduced dot (with one less electron) is renormalized 
and that new interaction terms appear beyond the conventional terms of the strong-coupling limit. 
The eigenenergies of this effective Hamiltonian are found to be in excellent agreement with exact 
numerical results of the original model in the limit of strong Kondo coupling. 
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I. INTRODUCTION 

The Kondo resonance, which emerges when a localized 
impurity spin interacts antiferromagnetically with a de- 
localized electron gas, has generated considerable inter- 
est over several decades by now.— 12, The observation that 
the Kondo resonance can be realized in the mesoscopic 
regime of quantum dots, in which many of the system 
parameters are experimentally tunable has led to much 
renewed interest over the last decade .-'^'^'^'^'^'i^iii'i^i'i^ 
This experimental work has been accompanied by much 
theoretical progress on the mesoscopic aspects of the 
Kondo problem^OeiiOSig^^^^^^^^ 

In the mesoscopic regime, the spin-1/2 Kondo impu- 
rity is typically represented by a small quantum dot with 
an odd number of electrons, while the delocalized elec- 
tron gas is realized by electrons in leads or in a large 
quantum dot. In this work we focus on the latter case, 
assuming a small and a large quantum dots that are cou- 
pled antiferromegnetically as in Fig.[l}i (see Ref.i for an 
experimental realization of such a setup). 

There are certain features that distinguishes the meso- 
scopic Kondo regime from the bulk limit. While the con- 
ventional Kondo theory assumes a continuum band of 
energy levels in the electron gas, the single-particle en- 
ergy levels in the large quantum dot are discrete. The 
discreteness and the dot-specific realization of these en- 
ergy levels become important when the Kondo tempera- 
ture Tk, the characteristic energy scale of the correlated 
Kondo resonance, is of the same order of magnitude or 
smaller than the average level spacing ,5, 17,19,23,24,27 
the conventional bulk Kondo model the electron-electron 
interactions in the electron gas are often neglected. How- 
ever, for the present double-dot setup, electron-electron 
interactions in the large dot can play an important role. 
In the following we assume the single-particle dynam- 
ics in the large quantum dot to be chaotic r^i^i^i^ in 
which case the dot is described by the so-called "uni- 
versal Hamiltonian" This Hamiltonian describes the 
low energy physics in a Thouless energy interval around 



the dot's Fermi energy. For a semiconductor quantum 
dot with a fixed number of electrons and in the limit of 
a large Thouless conductance, the electron-electron in- 
teraction is dominated by a ferromagnetic exchange in- 
teraction that is proportional to the square of the total 
dot spin. Despite its conceptual simplicity, this univer- 
sal Hamiltonian description was shown to yield a quanti- 
tative agreement's with experimental results measuring 
the statistics of the Coulomb blockade peak heights and 
spacingsj^ 

The effect of ferromagnetic exchange correlations on 
the Kondo resonance was first addressed analytically in 
the bulk limit^^ and, more recently, mean-field studies 
were carried out in the mesoscopic regime. a recent 
work, we studied this problem numerically and provided 
analytical results for the weak and strong Kondo coupling 
limits. Si We found that for weak Kondo coupling, the 
Kondo spin acts like an external magnetic field, assisting 
the ferromagnetic polarization of electrons in the large 
dot. In the case of strong Kondo coupling, the Kondo 
spin effectively removes one of the electrons of the large 
dot. We showed that this "reduced" dot with one less 
level and one less electron can again be described by a 
universal Hamiltonian but with a renormalized exchange 
constant. 

A central issue that was not addressed in our previous 
work concerns the nature of residual interactions in the 
reduced dot beyond the renormalization of its exchange 
interaction term. From the work of Nozieres'^ we know 
that a non- interacting electron gas turns into a Fermi- 
liquid when strongly coupled to a Kondo spin. The dom- 
inant effective interaction between the quasi-particles in 
this Fermi liquid is a repulsive interaction between spins 
of opposite orientation that are in close proximity to the 
Kondo spin. In the present case, the finite exchange in- 
teraction in the large dot leads to new effective interac- 
tion terms in the strong-coupling limit. To identify these 
new interaction terms, we follow a strategy that is sim- 
ilar to the one used by NoziereSf'^^ i.e., we perform an 
explicit perturbative expansion of the effective Hamilto- 
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nian of the reduced dot in the strong-coupHng Hmit. In 
the presence of exchange interaction, this strong-couphng 
expansion is significantly more involved. However, the re- 
sulting effective quasi-particle interaction contains only a 
few new terms. The analytical expressions we derived for 
these effective exchange-like interactions are validated by 
a comparison with a full numerical diagonalization of the 
original two coupled dot model. 

The outline of this paper is as follows: In section II 
we present the model of a spin-1/2 quantum dot that 
is Kondo-coupled to a large quantum dot (described by 
the universal Hamiltonian), and discuss the transforma- 
tion from the single-particle orbital basis of the large 
dot to a chain site basis, commonly employed in the 
strong-coupling limit. In Section III we discuss the limit 
of strong Kondo coupling and use a projection method 
to derive an effective Hamiltonian for the reduced large 
quantum dot with one less electron. In section IV we 
describe the evaluation of the eigenenergies of this effec- 
tive Hamiltonian, and in section V we compare the results 
derived from this effective Hamiltonian with an exact nu- 
merical solution of the original model. In section VI we 
conclude with a summary and discussion. 




2D 



•-■■■□^D □ Eh^A-D 

K 0^1 2 3 N-1 



Hko Hep 



FIG. 1: (Color online) (a) Schematic illustration of the system 
under consideration: a small quantum dot (left) with spin 
(Kondo spin) is coupled antiferromagnetically (coupling con- 
stant Jfc) to a large quantum dot (right) with spin S^. The 
large dot is described by the universal Hamiltonian, char- 
acterized by a ferromagnetic exchange interaction (coupling 
constant Js). The A*' single-particle energy levels in the large 
dot are distributed within a band of width 2D (half-filling). 
The average single-particle level spacing is fixed and given by 
S. (b) The large dot is represented in the site basis (squares) , 
in which Sk couples only to site 0. In the strong-coupling 
hmit, Jfc oo, it is useful to divide the Hamiltonian into 
three parts, Hko, Hep, H'd, according to the sites involved 
[see Eq. ©]. 



II. MODEL 

We consider a chaotic quantum dot that is coupled an- 
tiferromagnetically to a Kondo spin as realized, e.g., by 
a small quantum dot with an odd number of electrons. 
A schematic illustration of such a system is shown in 
Fig. [T^. In the following we will refer to the large quan- 
tum dot as the "dot" and to the small dot as the "Kondo 
spin." 



A. Hamiltonian 

In the limit of a large Thouless conductance, a quan- 
tum dot whose single-particle dynamics are chaotic is de- 
scribed by the universal Hamiltonian^^ 



JV-l 



H, = 



EE' 



(1) 



Here ajj ^ creates an electron with spin up/down [a — 
±1) in level n in an orbital single-particle level with en- 
ergy e°. We assume N spin-degenerate single-particle 
levels spanning a bandwidth of 2D = [N — 1) x J (J is 
the mean level spacing) . The second term on the r.h.s. of 
Eq. ([T]) represents a ferromagnetic exchange interaction 



(J, > 0) where Sd = i 



(x are Pauli 



matrices) is the total spin of the dot. In Eq. ([T]) we have 
ignored a constant charging energy term and a repulsive 
Cooper channel term. 

The dot is coupled antiferromagnetically to a Kondo 
spin {Sk - 1/2)^3 



H = Hd + JkSK- Sd(0) , 



(2) 



where Jk {Jk > 0) is the Kondo coupling constant and 
Sd{0) is the spin density of the dot at the tunneling po- 
sition r = 0. The dot spin density at position r is given 
by 



(3) 



where "01- (r) creates an electron with spin a at posi- 
tion r. In terms of the single-particle orbital wave 
functions 0„(r), the field operator is given by V-'J(r) = 
E^Jp^ 0„(r) ajj £r s-iid the local density of states of the 

dot is given hyil^ p{e) = En^o l<^40)P<5(e - 
with an average value of p « 1/{NS). 



B. Chain site basis 

The strong-coupling limit of the system in Fig. [1^ is 
more clearly described when the Hamiltonian in Eq. Q 
is rewritten in a different basis, known as the chain site 
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basis. This new basis is obtained by a unitary transfor- 
mation of the orbital basig- 

JV-l 
n=0 

such that site = corresponds to the tunneHng posi- 
tion r = 0, and the new one-body site Hamiltonian of the 
dot is tridiagonal, i.e., each site is coupled to its two near- 
est neighbors. Such a transformation is constructed by 
choosing c|j^q ^ = VjJ.(r = 0) and carrying out a Gram- 
Schmidt orthogonalization procedure.^ 

The chain site single-particle energies ej^ are given by 

the diagonal elements of Hq = J2n=oJ2a'^nai,a<,a 
when the latter is rewritten in the site basis. The off- 
diagonal matrix elements = and t* = t^^^-i de- 
scribe the hopping amplitudes between neighboring sites. 
A spin can be associated with each site, where the spin 
of site /i = is equal to the spin density at the tunnel- 
ing position, i.e., Sq = 8^(0). In the site basis, the Kondo 
spin couples only to a single site /i = 0. The Hamiltonian 
in Eq. ([2]) is now given by 

H = Ho-Js^l + Jk^K-so, (5) 

where the total spin of the dot is = 'l2^=o ^p- Here 
Hq is the one-body Hamiltonian of the dot in the site 
basis 

N-l 

with iJhop being the hopping Hamiltonian 

N-2 

Hhop = I] + h-C.) . (7) 

1-1 = (T 

The site basis formulation is particularly advantageous 
for the strong-coupling limit — > cx) when the site ^ = 
effectively decouples from the rest of the chain. Accord- 
ingly, we decompose the Hamiltonian in Eq. ^ into three 
terms (see Fig. [Ud for a schematic illustration) 

H = Hko + H'^i + Hep , (8) 

where Hko describes the Hamiltonian of the Kondo spin 
plus site n — 0, H'^ is the Hamiltonian of a "reduced" dot 
with — 1 sites /i = 1, . . . , A — 1 and Hep contains the 
remaining coupling terms. Writing = S'^ + Sq, where 
= J2p=i ^^'^ ^Pi'^ ^^"2 reduced dot, we have 

Hko = eg "0 - Js Sq + Jk • sq (9) 
K = H'o- JsSf (10) 
Hep = -2J,So-S;, + ^(io4,.Ci,,+ h.c.) . (11) 



Hq is the "bare" Hamiltonian of the reduced dot, given 
by expressions similar to Eqs. ^ and ([7]) but with the 
sums over fj. starting at ^ — I. 

Here and in the following, operators in the reduced dot 
space of A^ — 1 chain sites fi — 1, . . . , N — 1 are denoted 
by primed quantities. For such operators, the summation 
over sites /i starts from /i = 1 rather than from = 0. 



C. Site basis with good spin quantum numbers 

The Hamiltonian H in Eq. ([5]) is invariant under spin 
rotations and therefore conserves the total spin of the 
system (Kondo spin plus dot spin) Stot = Sk + ~ 
Sifo + S^. To take advantage of this symmetry, it is 
convenient to use a basis for which both the total spin 
5tot and the corresponding magnetic quantum number 
Aftot = 'S'tot good quantum numbers. 

There are different ways to construct a basis with 
good total spin, but one of them is particulary useful 
in the strong-coupling limit Jk ^ to,Js. To zeroth or- 
der in to/Jk and Js/Jk, we can ignore the couphng term 
Hep, in which case the subsystem of Kondo spin plus 
site decouples from the reduced dot. The Hamilto- 
nian Hko is easily diagonalized by coupling the spins 
Sk and Sq to S^o = + sq and by using Sk ■ §0 = 

If site /i = is singly occupied, i.e., uq — 1, this spin 
coupling will lead to either a singlet Sko = (lowest 
energy) or a triplet Sko — 1 (highest energy). However, if 
site /i = is empty (no — 0) or doubly occupied (tiq = 2), 
the spin at site and the corresponding Kondo coupling 
term vanish. This results in an unscreened Kondo spin 
in a doublet state {Sko = 1/2), the energy of which is 
intermediate between the singlet and triplet states. 

We now construct a basis of good total spin that 
also reflects the division into singlet/doulet/triplet man- 
ifolds. The eigenstates of Hko are characterized by 
SkojMko with Mko being the magnetic quantum num- 
ber of Sko ■ The eigenstates of the reduced dot Hamil- 
tonian H'j^ with N ~ no electrons are characterized by 
Ij'S'j^M^), where S'^, M'^ are the spin and spin projection, 
respectively, of the reduced dot and 7' denotes all other 
quantum numbers distinguishing between states of the 
same S'^ . We then couple the above eigenstates of Hko 
with the eigenstates of the reduced dot to form states 
with good total spin and spin projection quantum num- 
bers S'tot 7 Aftot • This basis of the coupled system is given 
by 1^0, Sko; l' , S'^; S'tot, Aftot)- To keep the notation sim- 
ple, we omitted the quantum numbers Sk — 1/2 and 
So = no(2 - no)/2. 

Spin selection rules determine the allowed values of the 
reduced dot spin S'^ for a given value of the total spin 
S'tot- We have S'^ = Stot for the singlet subspace, S'^ = 
S'tot ±1/2 for the doublet subspace, and S^ = Stot, S'tot ±1 
for the triplet subspace. 
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III. STRONG-COUPLING HAMILTONIAN 

The strong-couphng Umit is defined by ^ to- Since 
to ^ D ^ N6, this hmit corresponds to JkP ^ 1, where 
p — 1/{N6) is the average single-particle level density per 
site. In the strong-coupling limit, the lowest eigenstates 
of H are dominated by the singlet manifold. The bare 
singlet Hamiltonian (in the limit when Hep is ignored) 
is given by the reduced dot Hamiltonian H'^ with — 1 
electrons (except for a constant shift). However, virtual 
transitions between the singlet and doublet /triplet mani- 
folds add correction terms to this Hamiltonian. Our goal 
is to determine the resulting effective Hamiltonian for the 
reduced dot in the strong-coupling limit. 



A. Projection technique 

In the limit of strong but finite Kondo coupling, the 
above three manifolds (singlet, doublet and triplet) are 
coupled to each other. Specifically, the exchange cou- 
pling term in Hep couples the singlet and triplet man- 
ifolds, while the hopping term between sites and 1 
couples each of the singlet and triplet manifolds to the 
doublet manifold. To account for these couplings we de- 
fine projection operators Pg/Pd/Pt on the correspond- 
ing singlet /doublet /triplet subspaces {Pg + Pd + Pt = ^) 
and decompose the wave function ip = "05 + ^rf + V't 
accordingly. -^i The Schrodinger equation for the coupled 
system Hip — Eip can then be written as 



P=s,d,t 



(12) 



where each of the two indices a, j3 assumes any of three 
values {s, c?, i] and Haf) = Pa H Pp. 

In the strong-coupling limit, our system is described 
to zeroth order by the singlet Hamiltonian Hss, which 
contains the bare reduced dot Hamiltonian H'^ (for A^— 1 
electrons) and Hko (which assumes a constant value), 
completely decoupled from each other. Higher order cor- 
rections come from the coupling terms in Hep which lead 
to an effective "dressed" Hamiltonian of the reduced dot. 
This effective Hamiltonian H'^^ is formally determined 
by eliminating -tpd and ipt in Eqs. (|12p and by writing a 
single equation in the singlet manifold H'^^ips = Eips, 
where H'^^ is given by 

= Hss + Hst { E — Htt I Hts 



Hsd + Hst I E — Htt I Hi 



ltd 



X <^ E 



Hdd + Hdt yE — Httj Htd 
Hds + Hdt i^E — Httj Hi 



(13) 



The diagonal components Haa in the above equation 
are determined by evaluating Hko [Eq- ®] in each of the 
three subspaces a — {s, rf, t}. The coupling terms in Hc-p 
[Eq. (|lip ] do not contribute to these diagonal components 
with the exception of Sq-S^ which contributes to Htt only. 
We find 



Hs 



Ps 



3 Jfc 3 Js 



4 4 
Hdd = Pdie-'ono + H'^Pd, 



H'AP. 



(14) 
(15) 

Htt = Pt{el + ^-^-^+H',~ 2Js So • S',)Pt . (16) 

Contributions to off-diagonal components Hap with 
a ^ (3 originate in -ffcp- The hopping term in Hep changes 
the spin Sko by ±1/2 and can only couple the doublet 
to each of the singlet and triplet manifolds, while the ex- 
change term So ■ in Hep only couples the singlet and 
triplet manifolds. 



B. Expansion in the strong-coupling limit 

The effective Hamiltonian EP^^ and the construction of 
a good spin basis in the previous subsection are exact, in 
that no approximations were made beyond the original 
Hamiltonian H in Eq. However, the form of 

H°^ is not very useful for practical calculations. In the 
strong-coupling limit Jk ^ to ^ N_S. Since the exchange 
constant Js is typically below ~ 5, the condition Jk ^ 
Js is automatically satisfied in the strong-coupling limit. 
We can then expand H"^^ in the two small dimensionless 
parameters to/ Jk ^/{Jkp) and Js/ Jk- We will do so up 
to fourth order in these parameters, where the expansion 
terms are measured in units of Jk (this energy unit is set 
by the energy of the singlet). 

The starting point for this expansion is the unper- 
turbed singlet Hamiltonian Hss , the eigenbasis of which 
is given by \no = 1, Sko = 0; 7', S'^; Stot = S'^, Mtot)- The 
corresponding eigenvalues are 



'-{Jk + Js)+e'o + E'l"^ - Js Stot {Stot + 1) , (17) 



where Em'^ are the eigenvalues of Hq. These unperturbed 

eigenvalues, iJm' , are the limiting solutions to which the 
eigenvalues of the full Hamiltonian in Eq. ([13]) con- 
verge for Jfc — > cx). The differences between Em^ and 
at large but finite values of Jk are induced by the vir- 
tual transitions from the singlet to the doublet or triplet 
subspaces. These virtual excitations, in turn, give rise 
to effective interaction terms in the reduced dot, denoted 
by 5H°^ . The full effective Hamiltonian in the singlet 
manifod is then given by H'^^ = Hss + SH°^ . 

The effective interaction terms in SEP^^ must be con- 
sistent with charge and spin conservation.^ In particular. 
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SH'^^ must be a scalar operator in spin space (i.e., invari- 
ant under rotations in spin space) and invariant under 
time reversal. This restricts the possible terms that ap- 
pear in the effective Hamiltonian in the strong-coupling 
limit. 



Scalar one-body terms, i.e., hi and {J2cr ' 



,+ h.c.) 



lead to a renormalization of the one-body part of the 
reduced dot Hamiltonian H'^. Two-body scalars that 
can be constructed from the spin Si at site 1 and the 
total spin of the reduced dot are sf, Si ■ S'^ and 
S^^. The first sf = |ni(2 — hi) is the Nozieres term 
known from the conventional Kondo problem^^ (in the 
absence of exchange, Jg = 0) but the other two terms 
are new. The scalar triple product iS^ • (si x S^) (the 
imaginary i is necessary for time-reversal invariance) does 
not lead to a new term since iS'^ ■ (si x S^) ~ ~Si ■ S^, 
while a fourth order invariant is given by S^"*. Other in- 
variants such as S'^ ■ J2a- cr' (^1 o- ''"o-.tr' C2 „/l + h.c, S'^ 



h.c. and Si • J2a 



h.c. are allowed but, as we shall see, they 



cancel out in the effective Hamiltonian. 

We rewrite the effective Hamiltonian in (|13p as 
Hss + Y^t=i where 



5Hi 
5Hi 



Hsd 

Hst 
Hst 



[E - Hm) + Hdt Htd 
1 



H, 



ds ■) 



E-Htt 
1 



(18) 



H.. 



td 

E — Htt E — H, 



H,u + h.c. 



dd 



Hst — Htd — Hdt — Hts 

E — Htt E — Hdd E — Htt 



where 



A 
B 
C 



-H^^ + e'o (1 



n-Oj 



(20) 



— J,[-3/4 - 5tot(5tot + 1) + , (21) 
6Jk 



3jr 



7 Hdt H, 



dt ntd ■ 



(22) 



In Eq. p^ . we omitted the product terms AC, BC, C*^ 
since their contribution is higher than fourth order, while 
the contribution of AB + BA can be shown to vanish 
identically. 

To keep track of the various contributions for each of 
the SHi, we label them in the following by SHi.j. These 
terms are understood to act only in the space of the re- 
duced dot while the Kondo spin and the spin on site 
arc locked into a singlet. The operators SHij in the 
reduced dot space are obtained by taking a partial expec- 
tation value (. . .)s in the singlet state. The correspond- 
ing operators in the full space are given, respectively, by 
Ps{- ■ ■)sPs- III Appendix A we list several relations that 
are useful in calculating the expectation values of various 
operators in the singlet space. 

The most dominant contribution to Eq. (fTO|) arises 
from the unity operator term (in the round brackets). 
We find 



6 Hi, 



3Jk 



\HsdHds)s 



4 
3J/c 



hop hop 



(23) 



where we have substituted H^d by the hopping Hamilto- 
nian between sites and 1, 



H 



(0,1) 
hop 



h.c. 



(24) 



In the terms 6H3 and 5Hi above we have replaced Hdd + 
Hdt (e - Htt) Htd by Hdd (the neglected term gives 

contributions that are higher than fourth order in the 
expansion parameters). 

We next expand each 5Hi to fourth order in the pa- 
rameters t^/Jk and Js/Jk- Since the energy E is of 
the order J^, the fractions appearing in each term can 
be brought to a form 1/(1 — X) with X being small 
in the expansion parameters. We then approximate 
1/(1 — X) « 1 + X + X'^ . In the following we summarize 
the explicit calculation of each term. 



1. Evaluation of SHi 
For 6 Hi we find 

SHi ~ ^^Hsd {l + A + B + C + A^ + B^]Hds, 



and used Eq. (jA-7[) . Alternatively, HsdHds describes the 
spin transitions illustrated in Fig. [2]a,, and the result in 
Eq. P5|) can be derived by using Table U in Appendix A 
to sum up all the corresponding transition pathways. 

The term containing A in Eq. (fTO|) yields corrections 
that are second order in t^j Jk- Using the difference in 
the values of H'^^ ho in the singlet and doublet subspaces, 
and Eqs. we find 



6Hi,2 = -iT-r{HsdAHds)s 



4 



{e'o-4) + iei-£o)ni + 



^ H (*1 ^l'^^2a + 1^-' 



(25) 



(19) 



Except for an additional constant shift of 
~ (sTfe) l^ol^(^o ~ ^1)' ^ one-body operator 
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that can be incorporated into the unperturbed singlet 
basis by simply redefining the site energy and hopping 

amplitude ti in the unperturbed Hamiltonian Hq as 
follows^! 



The term involving C in Eq. (jl9p is given by 



IN 

3Jfc 



2 V 3Jfc 



{el - eg) , (26) 

2' 



(27) 



The term involving B in Eq. contributes only for 
a finite exchange interaction ( J^. ^ 0). We have 

(28) 

Using the identi ties (|X3|l - ((X^ . Eq. ^ can be simpli- 
fied to give Eq. (|A-15p in the Appendix. Using S'^ — Stot 
in the singlet subspace, we obtain 



SHi,3 = --^ {HsdBHds)s = 2 J, ( 



(29) 

We note that SHi g, is a spin invariant in the reduced dot 
space. 
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11,-1X1,-11 
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FIG. 2: Spin transition diagrams used to derive 
the efltective strong coupling Hamiltonian in Eq. (|53|l . 
All transitions connect two singlet states, character- 
ized by the quantum numbers Sko = 0, Sko.z = 0. 
The intermediate transition pathways involve combina- 
tions of doublet (Sifo = 1/2, Sko.z = ±1/2) and triplet 
{Sko = 1, Sko,z = ±1, 0) states. 



5Hi 



16 



{HsdHdtHtdHds 



(30) 



This term appears in the conventional Kondo problem 
(where Js = 0) and is known as the Nozieres term. 
Nozieres found'^^ that this term yields an effective in- 
teraction in the singlet-space that repels opposite spins 
on site 1. This term is induced by virtual transitions of 
the type singlet-doublet-triplet-doublet-singlet. Once 
we insert a triplet projection Pt in the r.h.s. of Eq. (|30p . 
i.e., we write the corresponding singlet expectation value 
as (HsdHdtPtHtdHds)s, we can replace both Hsd and Hdt 
by the hopping Hamiltonian H^^^'^^ between sites and 



hop 

1 [see Eq. mi Using Htl'^^PtKo '^ 



{HsdHdtHtdHds) s ~ 



(0,1) 



A-(0,l)rV(0,l) 

^hop -"hop 



/rr(0,l) rrk>',i; rrl>',-L; „ 
^-"hop -"hop -"hop 

v(o,i)ry(o,i 

'hop -"hop 



(0,1) ^(0,1), 
hop / 



KKo'^Htl^'^U'. (31) 



With the help of Eqs. jX^ . ((X^ . and dHU), we 

then find 



6Hi 



\to\ 



4 --2 



16 



|io|'''^i(2 ~ fii) 



(32) 

An alternative way to calculate (5-ffi,4 is to use the spin 
diagram in Fig. [2^. It can be reduced to the transition 
diagram in Fig. with the help of Table U in Appendix 
A. 

The Nozieres term ([5^ vanishes when site 1 is either 
empty (ni = 0) or doubly occupied {ni = 2) but is neg- 
ative for ni = 1, thus favoring a singly occupied site 1. 

The contribution from in Eq. (fT9|) is evaluated using 
Eq. (|A-13P and Eq. (IA-14P and leads to a constant shift 

^^1,5 = ^TV {Hsd-^Hds)s 

6Jk 



jj-J \h?[\h\'' + el+el-2s^e,] . (33) 



Finally, the contribution from in Eq. (fTO|) is found to 
be 



6Hi 



(. 



to\^J^{S'^ + 2s^-S'^) , (34) 



V3Jfc 

where we have used Eq. (jA-16p . 



2. Evaluation of SH2 

We next turn to the singlet-triplet transitions as de- 
scribed by 6H2 in Eq. (jlSp . The corresponding expres- 
sion for 5H2 is given by 

^ " (l + D + E + E^^Hts, (35) 



6H2 



Jk 
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where 



D 

E 



1 r 

ii 
Jk 



5'tot (5*101 



+ 2 So • 



(36) 
(37) 



In Eq. psp we omitted the terms and DE, which can 
be shown to vanish. 

The dominating term in Eq. (|35p is the one involving 
the unity operator. The corresponding term induces a 
spin transition as in Fig. ^Bp. Using Hst = —2JsPs{so ■ 
S'jPt and Eq. ((X^ we find 



6H2 



Jk 



{HstHts)s — — 



6'2 



(38) 



The same resuh can also be obtained with the help of 
Table HTl in Appendix A. The contribution SH2,2 induced 
by the term involving D in Eq. (|35[) can be simplified 
by using {H,tH'^Hu)s = J! KHqS'^- 
with the scalar operator Hq, we find 



Since commutes 



•Jk 



0. 



(39) 



The contribution SH2,3 induced by the term E in Eq. ((35|) 
can also be simplified since S^^ acts in the reduced dot 
space (and therefore has identical action in the singlet 
and triplet manifolds). The resulting expression gives 
rise to transition pathways as shown in Fig. [^t, the sum 
over which is further simplified using (|A-10|) to give 



SH2 



1 



HstEHt, 



(40) 



The last term 5H2a in Eq. ([55]) . containing E^ , is found 
to be (5iJ2,4 — — (16 Js/J|) (so • S^)'' and corresponds to 
the transition diagram in Fig. [5^. Since Hgt, as a scalar 
operator, commutes with S'^^, all other terms vanish iden- 
tically. Using Eqs. (|A-3|) this expression can be simplified 
to give 



5H2 



{Hst 



H, 



_ s 0/2 



(41) 



3. Evaluation of 5Hz 

Following Eq. the subsequent contribution, SH-^, 
involves transitions to both the doublet and the triplet 
subspaces 

5Hz « -^Hst{l + D + E)Htd[l + A + B)Hds+h.c., 

(42) 

where the operators A — E are the same as introduced 
above. Contributions involving products between {D, E} 



and {A, 13}, respectively, are higher than fourth order in 
the expansion parameters and therefore not considered 
here. 

The dominant contribution, 6H31 originates in the 
transitions shown in Fig. [2jl. Using Table [T] we can, 
however, simplify this transition diagram to the one of 
Fig. Os, for which we obtain 

SH3,i = ^{HstHtdHds)s + h.c. = ^liol'si • S;, . 

(43) 



To simplify the term SH3^2, involving D in Eq. (1421) . 
we make use of [Hst, Hq] = 0, leading to 
4 



5H 



3,2 



Hst D Htd Hds + h.c. 

oJi, 



(44) 



8Js 
34 



The diagonal matrix elements of SH3^2 (when evaluated 
in the eigenstates of the reduced dot Hq) vanish. Off- 
diagonal matrix elements enter in second order perturba- 
tion theory (for the effective Hamiltonian) and would lead 
to a correction ^ Js\to\*/Jk that is beyond the fourth 
order approximation. A convenient choice for is 
the average energy of the two reduced dot eigenstates 
(appearing in the corresponding matrix element). This 
choice leads to 

SHsa = . (45) 

The next term dHs^s, produced by E in Eq. ([i^ . 
can be brought into the form of Fig. [2}: (using Table [l| , 
where a sum of all occurring terms can be identified with 
Eq. (|A-4|) . yielding 

SH3,3 = A_{HstEHtdHds)s+h-c. 

^2^\to\H,-S'd+h.c. (46) 



^4 



The contribution (5i?3,4 induced by A in Eq. is 
calculated to be 



'5-ff3,4 = -Z-f2 {HstHtdAHds)s 



h.c. 



(47) 



'k 

32Js 



|topsi-S', - i?^) + h., 



Following the same arguments used for the evaluation of 
(5i?3_2, we find 6H3^4 = 0. 

The term B in Eq. (|42p gives rise to a transition dia- 
gram as in Fig. [2ji which can be reduced to a diagram 
as in Fig. with the help of Table IIIII Further simplifi- 
cations involving several of the equations in Appendix A 
yield 

<5iJ3,5 = A (Hst Htd B Hds)s + h.c. (48) 



3J,f 



Jk V3Jfe 



\tof{S'i + 2si-S'a)+h.c. 
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4- Evaluation of SH^ 



(to fourth order in t^j Ju and Jsj Jk) by 



To determine the terms contributing to (5iJ4 (up to 
fourth order), we can make the following approximations 
in Eq. (HID: {E-BuY^ « -1/Jfc and {E-HadY^ ^ 
— 4/(3Jfe). The resulting expression for (5i?4 corresponds 
to a transition diagram as in Fig. [2f. Using Eqs. (jA-ip . 
COj) . (IX^ and we find 



5i/4 K, — {Hst Htd Hdt Hts) 



(49) 



3 |toP S'i + 2si.S:, 



Alternatively, we can also obtain this results by using 
Table 1 to get a transition diagram as in Fig. which 
can be simplified by using Eq. (IA-4[) . 



+ Jk 



80 J,. 536 J2\ I top 



itor(3§i-i) 



4 



where is a constant 

, 3(Jfc + J,) 4\to\^ 
''° 4 3Jfe 



^) l^oP |tip + (e§)^ + (eD^ 



(53) 



(54) 



and i/^ describes a renormalized universal Hamiltonian 
for the reduced dot 



5. Additional terms 

In the above calculations, we have replaced the energy 
eigenvalue E in Eq. (fT5|) by its unperturbed values i?'*'^ 
as given by Eq. ([T7]). However, additional terms to the 
effective Hamiltonian are found when corrections to E^^"^ 
are included self-consistently. To the order we are in- 
terested in, it is sufficient to consider SHi^i and 5i?2,i 
[see Eqs. (^5)) and Eq. ([35]), respectively] as corrections 

to 



E ^ 



4|in 



3Jfc 



J2 

7" 'S'tot('S'tot 



1). 



(50) 



Higher order corrections to E in Eq. ((50)) lead to terms 
in the effective Hamiltonian that are higher than fourth 
order. 

Adding now the two correction terms from Eq. ([50)) to 
£;'(0) in Eq. ([20]) gives the following term 



\to\' 



4J, 



a/2 



(51) 



A similar correction to £"(°) in Eq. ([36]) yields: 



4 7^2 ^ t4 ^ 

rrr ^■^s I , |2 a'2 , s ^/i 



(52) 



Up to fourth order, corrections to E'^-^^^ do not lead to 
additional terms in 6H3 and 6H4. 



C. Effective Hamiltonian: a complete expression 



Collecting all the terms found in the previous section, 
the effective Hamihonian = Hss+Et=i IS given 



H'd 



•Js^d 



(55) 



Hq is a renormalized one-body Hamiltonian of the re- 
duced dot, obtained from the original one-body Hamil- 
tonian of the reduced dot, Hq, by redefining si and 
ti according to Eqs. ([26]) and ([27]) . respectively. This 
tridiagonal Hamiltonian can be rediagonalized Hq = 
J2nYicr ^ft'^Icr^no- to define new effective single-particle 
orbitals a|^^|0) and energies of the reduced dot. Jg is 
a renormalized exchange constant 



7 P P 

-/.^ 



Js |toP 

27 



(56) 



The most dominant contributions in (|56p are positive and 
thus lead to a stronger exchange interaction in the re- 
duced dot than in the original dot, Js > Js- Since the 
Kondo spin and the spin at site are coupled to a singlet, 
the spin of the reduced dot = Stot , and is thus con- 
served (i.e., S'^ = ^tot and M'^ = S'^_^ = Mtot are good 
quantum numbers). 

The effective Hamiltonian of the reduced dot contains 
several additional interaction terms, see Eq. ([55]) . The 
term proportional to (3s^ — 1) is the Nozieres term, 
known from the conventional Kondo problem (in the ab- 
sence of exchange, = 0)4^ The term proportional to 
Si-S^ is a new effective interaction in the reduced dot that 
is induced by the finite exchange interaction ( ^ 0) and 
describes an exchange interaction between the spin at site 
1 and the spin of the reduced dot. This exchange interac- 
tion is to leading order antiferromagnetic but depending 
on the particular values of Jk and Js it can also become 
ferromagnetic. The last term in ([53p is a four-body term 
but using S^'' = ^totJ it is easily evaluated in terms of the 
conserved total spin Stot- It can be combined with the 
renormalized exchange interaction in the reduced dot by 
defining an exchange coupling that depends on the total 
spin. 
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IV. EIGENVALUES OF THE EFFECTIVE 
HAMILTONIAN 

The effective Hamiltonian of the reduced dot is vahd 
to fourth order in t^/Jk and Js/Jk (when its terms are 
measured in units of Jk). To determine its eigenenergies 
to this fourth order, it is sufficient to solve H'^^ in first- 
order perturbation theory. Both H'^ and S^* are diagonal 
in a basis of good orbital occupations and good total spin 
of the reduced dot, while a second-order perturbation 
theory of the remaining interaction terms leads to terms 
that are higher than fourth order in the combined power 
of to/Jk and Jg/Jk- 

As required by first-order perturbation theory, we eval- 
uate the expectation value of H'^^ in the unperturbed 
basis, i.e., in the eigenbasis of the renormalized universal 
Hamiltonian iJ^ of the reduced dot. For simplicity, we 
will denote the good spin eigenstates of H'^^ by |^) and the 
corresponding expectation values by (. . Most terms 
contained in H'^^ [see Eq. ([55)) ] are diagonal in this basis, 
leaving only a few terms which require a special treat- 
ment. 

Consider first the evaluation of the Nozieres term. The 
calculation of (sf )j simplifies for the lowest eigenstate of 
H'^ at each given spin value 5^ = ^tot ■ Those eigenstates 
of H'j^ with = S'j^ have a maximal spin projection 
with only spin up electrons in singly occupied levels and 
thus have n-^ as good quantum numbers (in contrast to 
a general eigenstate of i?^, where only the orbital occu- 
pation numbers are well defined). For these states we 
have 

(T=± 
(T=± 

The occupation {fiia-}^ can be calculated from {fiia)^ — 

12^=1 l^inP ("-fio-)4; where U' is the unitary matrix (of 
order A'^ — 1) transforming between the renormalized 
single-particle orbitals of the reduced dot (with creation 
operators ot^) and the site basis states ^— 1, . . . , N — 1 

N-l 
n=l 

For the good spin eigenstates of H'^^ with S'^ = M^, the 
expectation value of Si • S'^; is given by 

(si-S:;)^ = (l + 5i)(si,,)e (59) 

C!=± h=l 




E I I I I I I I I I I I I I 3_ 

1 2 3 4 5 6 Jj^p 

FIG. 3: (Color online) Energy difTerence A_E between the es- 
timate based on the effective Hamiltonian (|53p and the exact 
numerical result for the ground-state energy in the subspace 
Slot = 3 (energy is measured in units of Jk)- The quantum 
dot contains 11 electrons in 11 spin-degenerate single-particle 
energy levels e° with an arbitrary random matrix realization 
of the single-particle spectrum, £°, but with non-fluctuating 
orbital wave functions 4>n{Q) ~ l/-/iV. The results, shown for 
an arbitrary but fixed value Js/5 = 0.52 (red symbols), be- 
have like ~ l/{JkpY (black solid line), expected for a strong- 
coupling expansion up to fourth order in l/{Jkp) in the limit 

JfeP -+ CO. 

Using Eqs. (|57|) and we can calculate the lowest 
many-body eigenenergy for each total spin value of the 
Kondo Hamiltonian ([2]) up to fourth order in to/Jk and 

Js/ Jk- 

V. COMPARISON WITH EXACT NUMERICAL 
DIAGONALIZATION 

To validate our expression ([55]) for the effective Hamil- 
tonian, we compare our analytical results for the many- 
body energies in the strong-coupling limit with an ex- 
act numerical diagonalization of the Hamiltonian ([2]) in 
a good spin basis scheme we developed previouslyj2I»22, 

As a first test, we compare results for the lowest energy 
of a given total spin (e.g., Stot — 3)- In Fig. [31 we show 
the difference AE (in units of Jk) between the energy 
determined from the effective Hamiltonian ([53[) and the 
energy found from exact numerical diagonalization of Q 
as a function of JkP (at an arbitrary, but fixed value 
Js/S — 0.52). This energy difference AE has to scale 
as oc l/(Jfcp)^, which is the next order of correction in 
l/{Jkp)" beyond the threshold of accuracy considered 
here. The results shown in Fig. [31 confirm this scaling 
behavior and thereby the accuracy and completeness of 
our effective Hamiltonian. Similar results (not shown 
here) are found for other values of Stot and for both even 
and odd number of electrons in the dot. 

It is interesting to study the ground-state value 
of the total spin S'tot- This quantity was studied 
theoretically^^i^i and can be probed experimentallyi^ 
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The ground-state spin S'tot undergoes successive transi- 
tions to higher values (known as the Stoner staircase) 
when the exchange coupUng constant Js is varied be- 
tween Js = and a value Jg S where the dot becomes 
fully polarized. The transition steps in the Stoner stair- 
case are shifted by the Kondo interaction. In Fig. [4] we 
show numerical results for the ground state spin diagram 
in the two-dimensional parameter space of JkP for 

a particular mesoscopic realization of the single-particle 
Hamiltonian of the dot. The exact spin transition curves 
(colored lines) that separate regions of fixed ground-state 
spin S'tot are monotonically decreasing for JkP ^ and 
monotonically increasing for JkP^ 1- Note also that for 
the particular mesoscopic realization chosen in Fig. lU 
some values of S'tot (e.g., S'tot = 1, 3) never become the 
ground-state values of the total spin in the weak-coupling 
limit. In contrast, the ground-state spin assumes these 
values in the strong-coupling limit. Our analytical re- 
sults in this limit, shown by the dashed lines in Fig. 21 
are in very good agreement with the exact numerical re- 
sults down to values of JkP ~ 2. The dotted black lines 
in Fig. 2] are the corresponding transition lines when we 
do not include interaction terms beyond the renormal- 
ized universal Hamiltonian of the reduced dot, i.e., when 
we assume the effective Hamiltonian to be just given as 
rj+H'^ [see Eq. (|53p ]. These dotted curves converge much 
slower to the full numerical solutions (colored lines) than 
the dashed curves determined from the effective Hamil- 
tonian ([55)1 . However, both the dashed and the dotted 
lines reproduce the monotonic increase of the exact tran- 
sition curves with Jk for JkP > 2. We conclude that 
this increase originates in the renormalization of the ef- 
fective exchange coupling constant in the reduced dot, 
Js ~ Js{l + Js/ Jk + ■ ■ ■), which is contained in the ap- 
proximations used for both dashed and the dotted lines. 
The renormalized exchange constant Js decreases with 
increasing Jk, which implies, in turn, that the spin tran- 
sition curves move upward with increasing Jk- 



VI. DISCUSSION AND CONCLUSION 

We have investigated the strong-coupling limit of the 
Kondo problem when the screening electrons reside in 
a large quantum dot that is described by the univer- 
sal Hamiltonian. The novel feature of this model, as 
compared with the conventional Kondo problem, is the 
inclusion of discrete level spacings and electron-electron 
interactions in the form of a ferromagnetic exchange in- 
teraction (that is part of the universal Hamiltonian). 

We have followed here a procedure that was originally 
proposed in Ref. [s^ for the conventional Kondo prob- 
lem in the absence of exchange correlations in the elec- 
tron gas. As pointed out there, one can find the effec- 
tive Hamiltonian at strong Kondo coupling (T <^ Tk) 
by considering the bare strong-coupling limit Jk ^ to- 
This bare strong-coupling limit is the one for which we 
have now provided a closed expression of all interaction 



J/5 

s 



0.9 



0.8 



0.7 



0.6 



T 1 — I I I I I I I 1 1 — I I I I I I I 1 1 — I I I I I I 

St t=5 

tot 



0.01 




J I 



FIG. 4; (Color online) Ground-state spin S'tot of the system in 
Fig. [T] with finite exchange constant Js and Kondo coupfing 
Jk- We consider 11 electrons in a dot with N = 11 single- 
particle levels, featuring an arbitrary random-matrix realiza- 
tion of the single-particle spectrum e°, but non-fluctuating 
orbital wave functions 4>n{0) = 1/y/N. Lines show the transi- 
tion curves separating regions of fixed Stot. The estimates 
based on our strong-coupling expansion [Eq. (|53p ] for the 
Hamiltonian (dashed black lines) are compared with full nu- 
merical solutions (colored solid lines). The dotted black lines 
are the transition curves obtained when the strong-coupling 
limit is described by a renormalized universal Hamiltonian 
but without additional interaction terms [i.e., H"^ = r/ + H'^ 
in Eq. (O]. 



terms up to fourth order in to/Jk and Js/Jk when the 
electron gas is described by the universal Hamiltonian. 
However, if the band width D of this electron gas is 
very large D ^ 5, the limit of strong coupling can be 
effectively reached at much smaller values of Jk than 
those of the bare limit. For such a system with large 
bandwidth, the strong-coupling limit corresponds to a 
Kondo temperature Tx that is larger than the system's 
temperature and average level spacing, Tk S> T, 6- An 
important insight in Kondo theory is that the effective 
Hamiltonians of both strong-coupling limits are related 
by a scaling analysis^ in which the Kondo Hamiltonian is 
renormalized by successive truncations of the band width 
D, leaving the low-energy physics unchanged. As the 
reduced band width D ^ {or equivalently T — > 0), 
the renormalized Kondo coupling constant Jk diverges.— 
The coupling constants of the various terms in the ef- 
fective strong-coupling Hamiltonian are typically deter- 
mined by fitting the lowest excitations of the effective 
Hamiltonian with those obtained by a numerical solu- 
tion of the full problem. For a detailed review of this 
procedure see, e.g., Refs. l^ lilHi^ . 
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By deriving in Eq. (|53p the effective Hamiltonian 77°^ 
for the bare strong-couphng Hmit, we have completed 
successfully the first step in our goal to understand the 
strong-coupling limit of the Kondo problem in the pres- 
ence of exchange correlations in the mesoscopic electron 
gas. We found that the exchange interaction in the uni- 
versal Hamiltonian gives rise to two new terms in H°^: 
a four-body contribution S^^ that can be absorbed into 
a spin-dependent exchange coupling in the reduced dot, 
and a new interaction term Si • that describes an ex- 
change interaction between an electron in the vicinity of 
the Kondo spin and the total spin of the reduced dot. 
This term is induced by the virtual polarization of the 
Kondo singlet involving excursions to both the doublet 
and triplet subspaces. Unlike the conventional Kondo 
problem [Jg — 0), this interaction is non-local as it in- 
volves the total spin of all sites of the electron gas (after 
the removal of an electron at site 0). 

It would be of interest to identify similar new inter- 
action terms in the low-temperature behavior of a cor- 
related Kondo state with a large band width D. Our 
numerical diagonalization method of Ref. [13 is limited 
to a rather small band width, e.g., = 11 levels for 
the results shown in Figs. [3] and [H For such small band 
widths, the bare and renormalized strong-coupling limits 
essentially coincide, and we could use our numerical diag- 
onalization method to validate the analytical derivations. 
For a numerical solution at larger bandwidths, a numer- 
ical renormalization group (NRG) technique^ might be 
useful. The challenge for NRG is the inclusion of non- 
local correlations induced by the exchange interaction in 
the universal Hamiltonian. With numerical solutions at 
hand it would be interesting to investigate whether the 
renormalization of the large band width Kondo problem 
induces any other "leading irrelevant" interaction terms 
around the strong coupling fixed point beyond those we 
have identified here. 
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APPENDIX A 

In this appendix we provide various expressions that 
are useful for deriving the effective Hamiltonian. The 
notation used here follows the convention of Ref. [isl 

The product of a hopping operator c\ c,co,cr be- 
tween sites and 1 and its hermitean conjugate can be 
expressed in terms of the occupation number (no, ni) 



and spin (so,Si) operators at these sites: 

4,<tCi,^4,^'Co,^' = ^o(l - ni/2) - 2so • Si , (A-1) 

^l,<T^O,aCl,a'^l,a' = - ^o/2) - 2So • §1 . (A-2) 



The spin raising and lowering operators, S± = Sx^iSy, 
satisfy, together with 5*^, the usual su(2) commutation 
relations 



[Sz,S±]—±S±, [S+,S-]^2Sz 



(A-3) 



while = Sl + ^{S+S- + S-S+). The su{2) commu- 
tation relations for the cartesian components of the spin, 

Si,Sj = i'^^.^ijkSk, can also be written in the form 

S X S = z S. The triple scalar product of two spin oper- 
ators, e.g., and Si is then given by 



(A-4) 



Other operator relations between and operators on 
site 1 are 



S'lc 



i± 



r 



/2 I ^'^1± 



Clzr'S^, ir T c^iS"^^^ , (A-5) 



d '^1± 



-t C;'2,iil± , it ql I st ql (A 



Useful relations involve the expectation values in the 
singlet state of observables at site 



{S0,i)s 
{S0,i ^Oj)s 



, 

1 „ 



= 1, (A-7) 

(A-8) 
(A-9) 

(A-10) 



where so,i is the i-th cartesian component of Sq and e^fe 
is the third rank antisymmetric tensor. 

We can also derive the following expressions for singlet 
expectation values of the form {Hsd ■ ■ ■ Hsd)s- 



{Hsd ho Hds)s 

{Had hi Hds)s 

{Hsd hi Hds)s 

{Hsd hi Hds)s 

and 

{Hsd Hds)s = 
{Hsd ^'d Hds)s 



= \to\' , 
= 2|topni 



(A-11) 
(A-12) 
(A-13) 

|to|'(l + ni ~2ni+ni_) , (A-14) 



- \to\'(s'j-Asi-S'dS 



I 0/2 

d 



(A-15) 



(A-16) 
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Matrix elements of various observables within and be- 
tween the singlet and triplet manifolds are listed in Tables 
I, II and III. 



{^llOllV'2) 


|0,0) 




|i,o) 




(o,o| 


1 


V2Si.- 


2 Si,, 


-^2 51,+ 


(1,-11 




1 + 2 Si,, 


-V2Si,+ 





(1,01 


2 Si,, 


-V2Si,- 


1 


-~V2Si,+ 




-V2Si.- 





-V2Si,- 


1 -2 Si,, 



TABLE I; Matrix elements {■ipi\Oi\ip2) of the operator Oi = 
-^hop ''-'^hop Vl^ol^- The corresponding states -f/)! {ip2) are hsted 
in the left column (top row), and are characterized by the 
quantum numbers Sko, Sko,,- 



TABLE 11: Matrix elements {ipi\so ■ S'j^\Tp2). Notation as in 
Table U 



(Vl|62|V2> 


10,0) 






(1,01 


— ^i-^'d^i- + ^i+^'d'^i+) 


(1,+11 





TABLE III: Matrix elements (^i|02|^2) of the operator 62 = 
Htl^^dHl^ol^ /\to\''. Notation as in TableU 



(V'llso ■ Sdl7/;2) 


10,0) 


11,-1) 


11,0) 


11, +1) 


(0,01 





-S'd.-/V2 


— 'S'd.z 


Sd,+/V2 


(1,-11 




~Sd,, 


S'd,+/V2 





(i,o| 




S'd.-/V2 





S'd,+/V2 


(1,+11 


S'd,-/V2 





S'd,-/V2 


+S'd,, 
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